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ABSTRACT 

We identify some of the most Hi massive and fastest rotating disk galaxies in the local 
universe with the aim of probing the processes that drive the formation of these extreme 
disk galaxies. By combining data from the Cosmic Flows project, which has consistently 
reanalyzed archival galaxy Hi profiles, and 3.6/im photometry obtained with the Spitzer Space 
Telescope, with which we can measure stellar mass, we use the baryonic Tully-Fisher (BTF) 
relationship to explore whether these massive galaxies are distinct. We discuss several results, 
but the most striking is the systematic offset of the Hl-massive sample above the BTF. These 
galaxies have both more gas and more stars in their disks than the typical disk galaxy of similar 
rotational velocity. The “condensed” baryon fraction, fc, the fraction of the baryons in a dark 
matter halo that settle either as cold gas or stars into the disk, is twice as high in the Hl- 
massive sample than typical, and almost reaches the universal baryon fraction in some cases, 
suggesting that the most extreme of these galaxies have little in the way of a hot baryonic 
component or cold baryons distributed well outside the disk. In contrast, the star formation 
efficiency, measured as the ratio of the mass in stars to that in both stars and gas, shows no 
difference between the Hi-massive sample and the typical disk galaxies. We conclude that the 
star formation efficiency is driven by an internal, self-regulating process, while fc is affected 
by external factors. Neither the morphology nor the star formation rate of these galaxies is 
primarily determined by either their dark or stellar mass. We also found that the most massive 
HI detected galaxies are located preferentially in filaments. We present the first evidence of 
an environmental effect on galaxy evolution using a dynamical definition of a filament. 
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1 INTRODUCTION 

How galaxies form and evolve remain open questions. Models re¬ 
produce global properties of galaxies well, partly by construction, 
but discriminative tests of plausible models are few and far be¬ 
tween. The greatest tests of theories and models often come from 
considering the extremes of parameter space. For example, some 
of the most demanding tests on the current paradigm of structure 
formation come from the lowest mass galaxies, which appear to 
be staggeringly underabundant (the “missing satellite problem”; 
[Moore et al.|1999[|Klypin et al.|1999| l. 

Here we focus on the opposite extreme by examining some of 
the most massive galaxies. Rather than considering their numbers 
as a population, we will examine their internal structure. While the 
least massive galaxies test certain aspects of the standard paradigm, 
namely the low mass end of the halo mass function, the nature of 
parent vs. satellite galaxy, and the role of reionization, giant galax- 
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ies help test the effects and results of a Hubble time of accretion, 
merging, collisions, and cannibalisms of dwarf companions (cf. 
IWang et al.|20T7> . 

In a theoretical framework of structure growth whose principal 
characteristic is hierarchical accretion, it is natural to suspect that 
the environment in which a galaxy resides plays a prominent role 
in its evolution. However, the success of models in which observ¬ 
able galaxy properties (luminosity, color, morphology) are based 
primarily on the galaxy’s stellar or dark matter mass (see for exam¬ 
ples the large literature on halo occupancy models, some examples 
includejPeacock & Smith|2060}|Berlind & Weinberg|2002[|Conroy| 
|& Wechsler|2009j l has been interpreted as supporting the view that 
environment plays a subdominant role in galaxy evolution ( |Kauff-| 
jmann et al. [2003[[Bundy et al.|2006| l. 

That is not to say that environment plays no role in such mod¬ 
els. The parent vs. satellite nature of a galaxy (cf. [Berlind et al.[ 
[2005|[Tmker et al.[2012[) is envisioned t o be key to explaining some 
critical observations ( [Peng et al.|2012[|Knobel et al.|2013[ l. How¬ 
ever, disentangling the physics involved is complicated greatly by 
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the different degrees to which environment correlates with various 
galaxy properties, such as morphology and color ( [Blanton et al.| 
|2003| ). And so, efforts to include environmental influences in mod¬ 
els increase the complexity and complicate interpretations of the 
data (cf. jWechsler et al.|2006[|Tinker et al.|2012[|Woo et al.|2013] l. 

Empirically, the role of environment in galaxy evolution is of¬ 
ten explored in the densest environments (see [Boselli & Gavazzi] 
|2006| for a review), which are necessarily the most complicated. 
Furthermore, the principal influence of the environment may not 
be direct, through phenomenon such as ram pressure stripping or 
tidal interactions, but rather indirect through either “assembly bias” 
( jSheth & Tor^|2004| [Gao et al.||2005| jWechsler et aTl|2006l l, 

which is purely a dark matter halo phenomenon in which close 
halo pairs tend to have earlier formation times than otherwise sim¬ 
ilar distant pairs, and “history bias”, in which the formation age 
of halos and their subsequent star formation history varies across 
environments ( jPe Lucia et aL|2012] l. Identifying clean tests of the 
influence of mass and environment on galaxy evolution is impera¬ 
tive. 


Defying the developed expectation that mass is universally the 
dominant determinant of a galaxy’s appearance and that massive 
galaxies should not contain cold gas and not be actively star form¬ 
ing (cf. jPekel & Birnboim|2006| jCattaneo et al.|2006| l is a class 
of massive, rotation dominated galaxy that includes HiZOA J0836- 
43, the most Hl-massive galaxy known (Sc, Mhi - 7.5 x 10^^ M©; 
[Haynes et al.|201l| t, the fastest known rotato r UGC 12591 (SO/a, 
Vmax = 506 km s“^; Giovanelli et al.||l986| l, and a few other gi¬ 
ant disks found in the recent 40% release of the ALFALFA survey 
( [Haynes et al.|201l] l. However, these are rare galaxies. Current de¬ 
terminations of the local galaxy Hl-mass function illustrate that sta¬ 
tistically, HiZOA J0836 -43 should not exist in the volume explored 
(probability of 3 x 10"^; Zwaan et al. 2005| l. Nevertheless, such ex¬ 


treme galaxies do exist, they are as massive as corresponding mas¬ 
sive early-type galaxies at low redshift ( [Bernard! et al.pOlT) , and 
they provide evidence by example against the hegemony of mass in 
determining a galaxy’s properties. 

Although such galaxies would clearly not be predicted in a 
model where mass alone determines a galaxy’s properties, we do 
not know whether more sophisticated models that are currently 
available successfully predict the properties of such extreme ob¬ 
jects. Because these objects are rare, global statistical tests of the 
models (cf. [Obreschkow et al.[2013[ l may not recognize the absence 
of such galaxies. It is therefore critical to complement large-scale 
comparisons, which offer the powerful advantage of statistics, with 
carefully crafted samples that aim to challenge the models at the 
margins. 

These unusual galaxies must have experienced a strikingly dif¬ 
ferent evolutionary history than that of massive early type galaxies, 
whether that means fewer or no large accretion events, a different 
angular momentum accretion history, and/or a distinctly different 
star formation history. They are likely to be the nearest analogs of 
the “classical” view of galaxy formation ( [Eggen et al.[1962[ ), where 
galaxies form smoothly within an overdensity. In fact, to produce 
disk-dominated galaxies numerical simulations of spiral galaxies 
set in a cosmological context have had to select to model systems 
in underdense environments, a reflection on what is now under¬ 
stood to relate to the details of the angular momentum accretion 
history (cf. [Governato et al.[2007[[Sales et al.[2012[ >. Therefore, the 
importance of environment in the genesis of these galaxies is hard¬ 
wired into detailed physics-oriented models, in opposition to what 
is often done in statistically-oriented theoretical treatments. 

Complementing the morphological differences between early 


and late type giant galaxies, there are differences in star formation 
histories. Late type giant galaxies presumably provide a measure 
of the quiescent mode of star formation and a contrast to what is 
likely to have occurred in early type systems. Empirically, we need 
to determine whether these differences are reflected in the overall 
star formation efficiency, which can be defined either relative to the 
total “condensed” baryonic content (the mass in stars and cold gas 
within the disk) or to the total mass. The former measures the frac¬ 
tion of all “condensed” baryons that formed stars, while the latter 
measures the fraction of the total mass converted to stars and so 
includes any variations in the “condensed” baryon fraction among 
galaxies. 

At this point, we pause to consider the ill-defined, and of¬ 
ten ill-measured, description of a galaxy’s mass. Theoretically, the 
preference would be to describe galaxies by their total mass, where 
total refers to both dark and baryonic mass, or halo mass, where 
halo refers to the dominant dark matter component. Even there, 
these are moving targets because of the definition of an outer ra¬ 
dius, such as a virial radius. However, such definitions are com¬ 
pletely impractical in an empirical sense because we rarely have 
any direct measurements of an individuals galaxy’s total or halo 
mass. Instead, investigators either default to estimating the bary¬ 
onic mass, which can be problematic in its own ways (as discussed 
next), or to assuming that the ranking between baryonic (or, typ¬ 
ically, more accurately stellar) mass and halo mass is nearly 1:1 
(see, for examples, halo abundance matching schemes; [Berlind &| 
[Weinberg|2002[ l. The latter approach has worked well in reproduc¬ 
ing statistical measurements of galaxy populations from theoretical 
dark matter distribution models, lending support to the hypothesis 
that the rankings are not grossly different than 1:1. We will there¬ 
fore measure baryonic masses, and presume that these are not only 
measurements of the baryons in these systems but also, in a ranked 
sense, measurements of the total mass. 

The baryonic mass measurements come with their caveats as 
well. It is usually the case that only stellar mass, through luminos¬ 
ity and colors, is estimated and corrections for baryons in other 
phases are either ignored or treated in the mean. We will describe 
a somewhat more accurate process in which we use IR magnitudes 
that have been carefully calibrated to stellar masses and comple¬ 
ment those with measurements of the gaseous baryons based on Hi 
observations. This approach still ignores baryons in hot phases and 
there is a well-known baryon shortfall ( [Bregman][2007[ ), so these 
masses are still not the full picture. However, to be specific, we will 
measure the stellar plus gaseous masses, and use these liberally as 
a proxy for the total galaxy mass. This assumption is broadly valid, 
not worse than a factor of two in total mass, as evidenced by the 
scatter about an apparently constant baryon fraction resulting from 
studies of the baryonic Tully-Fisher relation ( [Zaritsky et al.[2014b] l 
for galaxies spanning many decades in baryon mass. 

The Cosmic Flows project, which is consistently remeasur¬ 
ing all available Hi line widths, provides a uniform, unprecedented 
dataset from which to draw these massive rotating and Hl-massive 
galaxies. In combination with the infrared imaging provided by the 
Spitzer Space Telescope and the calibrated conversion between 3.6 
photometry and stellar mass ( [Eskew et al.[[2012] |, we are now in 
a position to compare gaseous, stellar, and total masses for these 
rare rotationally supported giant galaxies. In §2 we present the data 
and present our two subsamples of extreme galaxies in §3. In §4 
we examine the properties of these galaxies relative to the bary¬ 
onic Tully-Fisher relation and draw inferences regarding the “con¬ 
densed” baryon fractions and star formation efficiencies. We sum- 
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marize our conclusions in §5. Where needed, we adopt Hq = 70 km 
s"^ Mpc"^ 

2 THE DATA 

2.1 The Cosmic Flows Program 

Since 2009, the Cosmic Flows project (CF) has gathered all the dig¬ 
ital Hi spectra available from public archives of the largest radio- 
telescopes worldwide and measured them in a consistent way. Two 
sub-projects of CF, at Green Bank in the USA and at Parkes in Aus¬ 
tralia ( [Courtois et al.| 20 lT] ), complete the archives for targets with¬ 
out previous observations adequate for the Tully-Fisher method of 
distance measurements ( |Tully & Fisher] 1977] l. The latest update of 
this catalog is published in [Courtois & Tull^ ( |2015| l. 

The main goal of CF is to map the all-sky peculiar velocity 
field locally and to uncover, at high spatial resolution, the underly¬ 
ing dark matter distribution. In that cause, they have measured tens 
of thousands of galaxy line widths with a new method described 
by [Courtois et al.| ( |2009| ) and [Courtois et al.[ ( [2011| ). Briefiy, for 
each galaxy, they evaluate the line width parameter Wm50, which 
measures the Hi profile width at 50% of the mean fiux within the 
velocity range encompassing 90% of the total Hi flux. That quantity 
is transformed into the more physically motivated parameter 
which is the same width corrected for the slight relativistic broad¬ 
ening and for broadening due to finite spectral resolution, corrected 
for inclination and so is effectively twice the maximum rotation ve¬ 
locity. Technical details of this process are provided by [Courtois 
et al.[([2009[) and [Courtois et al.[d2011[), and reviewed in [Tully & 
Courtois|j2012t . - - - 

The result is a catalog of Hi measurements of unprecedented 
size (14221 galaxies) and consistency. They use the ratio of the sig¬ 
nal level at 50% of the mean flux to the noise measured beyond the 
frequency extremities of the signal to parametrize the line width 
uncertainty. Only profiles with uncertainty estimates smaller than 
or equal to 20 km s"^ are retained (more than 75% of the 14221 
galaxies, exactly 10733 galaxies satisfy this criteria) after a supple¬ 
mentary visual inspection. This catalog is available for public use 
at the Extragalactic Distance Database (FDD) websit^ where sev¬ 
eral other parameters, such as the integrated Hi line fluxes for those 
lines with a flux calibration error better than about 10 to 15% and 
the average heliocentric velocities, are also available and described 
( [Tully et al.|2009[ l. In addition to this “all-digital” catalog, we also 
use the “pre-digital” catalog of measurements at 20% of the Hi line 
peaks. The relation 

\% 5 o(all - digital) = (1.015 x W 2 o(pre - digital) - 11.25) (1) 

connects the two catalogs { [Courtois et al.|200^ . In total, we have 
16,124 galaxies with coherent Hi measurements and 12,189 with 
sufficiently good quality for distance measurements. 

2.2 The Infrared Data 

We use 3.6yL/m magnitudes measured from data obtained using the 
Spitzer Space Telescope ( [Werner et al.[2004| | and the IRAC camera 
( [Fazio et al.|2004[), from a number of different observing programs 
(such as, [Dale et al.[2009[[Sheth et al.[2010|[Sorce et al.[201^ . For 

Tully-Fisher studies, within the CF sample, there are 6,007 galax¬ 
ies with sufficiently good Hi linewidth measurements, crv^mso < 20 

^ http://edd.ifa.hawaii.edu; catalog “All Digital HI” 


km s"^ and inclinations greater than 45°. Among these, 2,493 have 
been observed with Spitzer. We ado pt total magnitude s measured 
either us ing the Archangel pipeline (|Schombert|2007[ ) or the 
pipeline ( [Munoz-Mateos & et al.|2015| ). Distances, needed to cal¬ 
culated absolute magnitudes, are presumed from the Hubble flow 
using Ho = 70 km s“^ Mpc“^ and the recessional velocities from 
the Hi catalog. 

We use the 3.6pm photometry because there now exist sim¬ 
ple, and robust, conversions to stellar mass ( [Eskew et al. [[20121 
[Meidt et al.[2014[ l that have been vetted with comparisons to SDSS 
SED-derived stellar masses ( [Cybulski et al.[|MT4] ). [Eskew et al.[ 
\2012\ give the conversion between stellar mass and infrared flux 
as M* = 10 ^-^^F 3 . 6 (q^)^, where D is the distance to the galaxy in 
Mpc, F 3.6 is the flux in Jy, and M* is expressed in solar masses. 
A slightly more precise estimate can be obtained by including the 
4.5pm flux, but both [Eskew et~^ ( [2012[ ) and [Querejeta & et al.[ 
( [2015[ ) find only a weak dependence between and infrared color. 
The largest uncertainty in the derived stellar mass remains the over¬ 
all normalization due to uncertainties in the adoption of a specific 
stellar initial mass function (the [Eskew et al.[ ( [2012[ ) calculation 
adopts a Salpeter IME). 


3 GIANT DISK GALAXIES 

To define our galaxy subsamples, we search both for those galaxies 
with the largest Hi masses and for those with the largest rotational 
velocities. Due to our use of the Hi database, we may miss gas- 
poor versions of the latter, although such galaxies are likely to be 
early type and therefore not the pristine test cases we are searching 
for. Because of broad correlation between gas mass and total mass, 
these samples have some objects in common. 

3.1 The most Hi massive galaxies 

We compute Hi masses (Mhi) from the publicly available integral 
fluxes and average heliocentric velocities ( [Courtois et al.[200^ us¬ 
ing 

Mhi = 2.36 X lO^D^FI (2) 

where D is the distance in Mpc and FI is the integrated flux in the 
Hi line in Jy km s“^. 

The distribution in log Mhi for the 10733 galaxies for which 
we can calculate the necessary quantities covers the range of 5.53 to 
10.72, with a median of 9.62, corresponding to a range of 3.4 x 10^ 
to 5.2x lO^^M©. PGC8681/UGC01752 is the most Hi massive (Fig¬ 
ure!^ galaxy in the full sample, but it is not in our subsample be¬ 
cause we lack Spitzer photometry for it. PGC17625, with only a 
slightly lower Hi mass (4.5 x lO^^M©), tops our list. In Figurej^we 
show the distribution of Mhi for the 100 most Hi massive galax¬ 
ies that also have existing Spitzer photometry. Choosing the top 
100 , which was set only because it is a simple round number, cor¬ 
responds roughly to selecting the top 5% of the available sample. 
One galaxy (PGC90167) was removed from our subsample b ecause 
even though it was observed with Spitzer, the S"^G survey ( [Sheth 
[et al.|2010| [Munoz-Mateos & et al.|2015] l was unable to provide a 
reliable magnitude measurement. 

We distinguish this sample from an Hi rich sample, where in 
that case one would presumably select on either the ratio of gas to 
stellar mass, Mhi/M^, or to dark matter mass, Mhi/Mdm- Although 
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Figure 1. HI profile of PGC0008681 as observed at Arecibo, and B band 
image from DSS2 survey at the size of Arecibo beam (3 arcmin). This is 
the largest HI mass galaxy recorded in the full available sample, but is not 
in our subsample due to lack of Spitzer photometry. 


such a sample is also of interest, that selection depends on the mod¬ 
eling of M* or Mdm and therefore complicates any interpretation. 


3.2 The fastest rotators 

Our rotational velocities come from the W^x parameter: 


^ W^J2 
sin(/) 


(3) 


where Vrot is the galaxy rotational velocity and i the inclination to 
the line of sight. 

Among the full sample of galaxies (not yet implementing the 
Spitzer imaging requirement), the rotational velocities range from 
11 to 514 km s-‘ for PGC71392AJGC12591 (see Figurewith a 
median of 145 km s"'. Again the top ranked galaxy is not in our 
final sample due to the lack of infrared photometry, and so the top 
ranked galaxy in our sample is PGC17625 with a rotation veloc¬ 
ity of 450 km s“^ The distribution of rotational velocities of our 
sample is also illustrated in Figure 


4 DISCUSSION 

4.1 The baryonic Tully-Fisher relationship 

Because in any sample there will always be some galaxies with 
the largest Hi masses or rotation velocities, we require a fiducial 
against which to compare these galaxies to the average galaxies. 


Figure 2. Distribution of galaxies in our subsamples selected by rotational 
velocity and HI mass. We highlight the HI mass selected sample using open 
blue circles and the rotational velocity selected sample using filled red cir¬ 
cles. We will continue this color scheme throughout the paper to differenti¬ 
ate the samples. There are 26 galaxies in common between the two samples. 


Are they simply scaled up versions of lesser galaxies or are they 
otherwise distinct? 

This class of question can be addressed using galaxy scaling 
relations. For rotationally supported galaxies the most commonly 
used such relation is that identified by |Tully & Fisher| ( ri977| l, which 
relates rotational velocity to optical luminosity. One can then use 
such a relation to ask such questions as whether these galaxies have 
unusual luminosites for their rotational velocities. 

Since the original work, this scaling relation has been ex¬ 
tended in one key way that has been shown to be particularly nec¬ 
essary for lower mass galaxies, namely the substitution of the bary¬ 
onic mass in place of the luminosity ( Freeman] 1999 [Walker 1 1 999] 


[McGaugh et al.|2000||Verheijen|200l[|Geha et al.|2006| l. For our 

purposes, we wil l use the baryonic Tu lly-Fisher (BTF) measured 
for S"^G galaxies ( [Zaritsky et al.|2014b| ) as our fiducial because the 
data used to establish that relation are most comparable to the data 
presented here. 

To implement the BTF we require measurements of the stel¬ 
lar and gaseous masses. We described in §2.2 how we obtain M*. 
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Figure 3. HI profile of PGC0071392 as observed at Arecibo, and B band 
image from DSS2 at the size of Arecibo beam (3 arcmin). This is the largest 
rotational velocity recorded in the full available sample, but is not in our 
subsample due to lack of Spitzer photometry. 

Obtaining a measurement of the gaseous mass, Mqas. involves cor¬ 
recting the Hi mass for the missing He, other “metals”, and molec¬ 
ular gas. This is done here as described by |Zaritsky et al.|f2014b| ) 
to enable a direct comparison to their results. The baryonic mass is 
then simply M* -h Mqas- As noted by those authors, this quantity is 
more accurately referred to as the “condensed” baryonic mass be¬ 
cause it only represents the baryons that have settled into the disk 
of the galaxy and is missing the possibly substantial baryonic com¬ 
ponents in the halos of these galaxies. 

In Figure|^we present our first comparison of the properties of 
the giant galaxies to the more general population. As shown, both 
the fastest rotators and most Hi massive galaxies lie generally on 
the BTF relation, although some differences exist. First, although 
slightly difficult to notice from this Figure but easier to spot in Fig- 
urej^where we plot the residuals from the mean BTF (ATFBaryons)» 
the fastest rotators fall somewhat below the mean trend as depicted 
by the red line in Figure Second, and much clearer, the Hl- 
massive sample lies systematically above the line. The one striking 
outlier, PGC90167, has a T-Type = -2. Because it is such an early- 
type galaxy, its dynamical support is probably not dominated by 
rotational support and therefore the galaxy should not be expected 
to satisfy the BTF (it has a low \rot of 45 km s"^). We remove this 
one galaxy from further consideration, leaving both the fast rotator 
and Hl-massive samples with 99 galaxies, and proceed to discuss 


Figure 4. Baryonic Tully Fisher. Comparison of the BTF for the S^G sam¬ 
ple and our subsamples of extreme disks. The red filled circles represent the 
fast rotators, while the blue open circles represent the Hl-massive galaxies. 
The red solid line is described in detail by |Zaritsky et al.| ( [2()14b) but in brief 
represents a model where a fixed fraction (~ 0.4) of the baryons in a halo 
of circular velocity of \rot, which has the cosmological baryon fraction, are 
settled in the disk. As shown by those authors, this simple model is consis¬ 
tent with the empirical best fit line. The dotted lines represent factors of two 
change in that fixed fraction. 


the two subsamples and their location relative to the mean BTF re¬ 
lation. 

There are a variety of potential explanations for why the fast 
rotators fall below the mean BTF: 

First, this result could be a manifestation of a Malmquist-like bias. 
Because we select this subsample using Vrot, errors that artificially 
inflate Wrot will preferentially contaminate our sample with galax¬ 
ies that appear to be “baryon poor” for their rotation velocity. The 
mean of the distribution (Figurej^ is displaced by ~ 0.15 dex, sug¬ 
gesting that we would need to have a similar magnitude error in Vrot 
for this to be a plausible explanation. That level of error translates 
to an error of ~ 40% or over 100 km s"^ for Vrot = 300 km s"^ 
Given our demands on precisely determined Hi widths and high 
inclinations, this level of error seems unlikely. 

Second, inclination errors, in the form of an underestimation of the 
true inclination, could lead us to preferentially place objects to the 
right of the mean BTF. We test for this effect by searching for a 
correlation between inclination and ATFearyons- There is only a very 
weak positive correlation (positive is in the necessary sense) that 
is entirely consistent with arising at random (can be rejected as 
arising randomly with only 37% confidence). We conclude there 
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shown in Figure]^ for galaxies with Nrot > 2. Upper panel shows t he resid¬ 
ual distribution for the S^G sample from |Zaritsky et aLH2014b) , middle 
panel for our subsample of fast rotators, and lower panel for our subsample 
of Hl-massive galaxies. Also shown in the lower panel, with the opposite 
hashing is the subpopulation of Hl-massive galaxies that are also in the fast- 
rotator sample. The small horizontal bars placed above each histogram rep¬ 
resent the mean of each distribution and the length represents the standard 
deviation of that mean. 


is no evidence that inclination is systematically affecting the BTF 
residuals of the high Wrot galaxies. 

Third, at these high Wrot values we may have a larger fraction of 
galaxies for which pressure supported components play a larger 
dynamical role. |Noordermeer & Verheijen| p007| ), when examin¬ 
ing the high mass end of the Tully-Fisher (TF) relation, found that 
a number of their galaxies were SO/Sa types (as well as finding a 
similar displacement of these massive galaxies below the general 
TF). When dealing with such galaxies a number of investigators 
( [Burstein et al.| 19^ [Weiner et al.|2006l|Kassin et al.|2007[|Zarit-| 

|sky et al.|2008| ) have suggested combining the use of Wrot and veloc¬ 
ity dispersion, cr, into a single kinematic term that measures the full 
dynamical support. While this may be appropriate when observing 
stellar components, it is not for gaseous measurements (there is 
negligible pressure support of the gas component and, even if there 
were, such a correction would move these galaxies even further to 
the right in the Figure). If one does attribute the dip below the BTF 
for this subsample to the prevalence of early-type galaxies at these 
rotation velocities, the cause of the dip would have to be a general 
failure of (or departure from) the linear BTF, rather than from the 
omission of cr in the evaluation of the dynamical support. 

Fourth, the intrinsic BTF may deviate from a straight line. Numeri¬ 


cal simulations ( |Aumer & White|2013| ) also find that the most mas¬ 
sive galaxies tail off below the BTF ( [Zaritsky et al.||201^ and 
so this phenomenon may be a natural result, for which an intu¬ 
itive explanation might provide significant insight to the nature of 
galaxy formation. As a potentially interesting aside, massive early 
type galaxies also appear to fall below their analogous scaling re¬ 
lation, the Fundamental Plane ( [Bernardi et al.|[20TT| l. As for the 
early types, for which curvature in the scaling relation is established 
( [Zaritsky et al.|2006[ [Bernardi et al.|2011| ), a unified scaling rela- 
tion ( [Zaritsky Qi al.[|2008[ ^ritskypoT^ would imply the same 
phenomenon for late types. 

Finally, the position of these galaxies may indicate that the stellar 
mass has been underestimated in these systems relative to that of 
the lower \rot systems. This error could arise because the stellar 
initial mass function (IMF) is more bottom heavy in the faster ro¬ 
tators. Although this hypothesis may seem farfetched, there is em¬ 
pirical evidence for this type of behavior among early type galaxies 
( [van Dokkum & Conroy|2010| [Cappellari et al.|20T^ , where the 
IMF becomes progressively more bottom heavy as one considers 
galaxies with larger velocity dispersions. Evidence for IMF vari¬ 
ations even extends to Local Group clusters ( [Zaritsky et al.|201^ 
[2014a[ ), and so the possibility that such variations exist among disk 
galaxies as well should not be quickly dismissed. 

As already mentioned, the more striking departure from the 
BTF is seen for the Hl-massive subsample. Here the departure, 
shown clearly in Figures and can be interpreted to mean that 
these galaxies have a larger fraction of their cosmologically ap¬ 
portionment of baryons within their disks, as stars or gas, than the 
typical disk galaxy. As with the fast rotator sample, there are a set 
of potential systematic errors that could also produce this result: 

First, one might wonder whether this could be the result of a 
Malmquist-like bias. We select systems with large HI masses, and 
therefore errors in the measured HI mass could both help popu¬ 
late the sample and make it appear as if these galaxies are Hi rich 
for their rotation velocity. However, these galaxies also have larger 
stellar masses than typical, as shown in Figure where we com¬ 
pare to a stellar mass only version of the BTF (renormalized to 
produce a mean departure of zero for galaxies with log Wrot > 2 
in S"^G). Because the stellar masses, and this stellar version of the 
BTF, are entirely independent of Mhi, we conclude that the depar¬ 
ture of these galaxies from the BTF is not a result of errors in the 
Hi measurements themselves. 

Second, errors in the adopted distance would affect both the 
gaseous and stellar masses because distance plays a role in convert¬ 
ing between observables and masses. The mean shift between either 
the gaseous or stellar masses and the fiducial of the S"^G sample is 
~ 0.22 dex (see Figuresand [^. Because both mass estimates are 
proportional to this offset corresponds to a systematic offset of 
0.11 dex in D or alternatively a 30% error. We know such an er¬ 
ror cannot arise from Hq given the current degree of uncertainty in 
the Hubble parameter ( [Freedman & Mador^[2010[ l. Peculiar mo¬ 
tions could be as large as several thousand km s"^ in th e richest 
galaxy clusters, but are typically a few hundred km s"^ ([Peebles 
[T?76l [K aiser n^ , particularly in the low density environments 
that these galaxies reside in. One might dismiss peculiar motions 
on the basis that these should result in increased scatter rather than 
an offset, but recall that we select the most Hi massive galaxies, 
and therefore are predisposed toward galaxies with overestimated 
distances. However, even with peculiar motions that are as large as 
500 km s"^ the induced error is < 30% in 97 out of the 99 galaxies 
in the Hl-massive sample. We conclude that distance errors are not 
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responsible for the systematic offset of Hi massive galaxies from 
the BTR 

Third, inclination errors will affect \rot- However, unlike in the pre¬ 
vious discussion, inclination does not play a role in the selection 
of these galaxies (other than the requirement that they be inclined 
by more 45°). As before, we test for the influence of inclination by 
searching for a correlation between measured inclination and the 
offset from the BTR Again we And no significant correlation (the 
probability of the observed correlation arising at random is 18%). 
We conclude that inclination errors are not driving the BTR residu¬ 
als. 

Rinally, the adopted relation between n rot and halo mass may be in¬ 
correct. We have the taken a mean trend from |Bullock et al.| ( |2001| ), 
normalized using the Milky Way (see [Zaritsky et al.||2014b| ), but 
different types of galaxies will have different degrees of adiabatic 
contraction, and therefore different relations between the inner ro¬ 
tation curve, Wrot^ and the characteristic halo circular velocity, Vc. A 
factor of two overestimation in the fraction of baryons within the 
halo that have settled into a disk, /c, would arise if we had under¬ 
estimated the halo mass of galaxies in our subsample by a factor of 
two, or alternatively underestimated v^ by Because it is diffi¬ 
cult to imagine a scenario by which Wrot is smaller than that due to 
the dark matter, we would posit that the problem, if it exists at all, 
is that the v^ of the typical galaxies are overestimated by a factor 
of V 2 . Ror a typical galaxy in our sample with Vrot ~ 180 km s"^ 
the necessary error would require the halo to have Vrot ~ 127. Such 
a steep drop between Wrot and v^ is not found in dynamical analy- 
ses ([Zaritsky & White] 1994| ) or gravitational lensing ( [Rischer et al.| 
|2000| ). While the conversion between Wrot and v^ has its associated 
systematic uncertainties, we conclude that it is unlikely to globally 
cause a factor of two increase in the apparent fc of the Hi massive 
sample. 

In the upper panel of Rigure|^ we calculate how the BTR off¬ 
set translates into differences in the level of baryon “condensation” 
in galaxies (see |Zaritsky et al.|2014b| for the original discussion of 
this term and the calculation relative to the BTR) as a function of 
Wrot- In other words, the Rigure shows the fraction of baryons — 
stars and gas — apportioned to a given halo at the level of the uni¬ 
versal baryon fraction that settle into the disks of these galaxies. In 
the lower panel of the same Rigure, we plot the star formation effi¬ 
ciency, or alternatively the fraction of the “condensed” mass that is 
turned into stars, = M=,/(M=, -h Mgas)- 

The Hi massive galaxies are clearly distinct from the typical 
disk galaxy in that the fraction of condensed baryons, which in 
some cases reaches nearly 100% of the cosmological baryon frac¬ 
tion (meaning that all of the baryons expected in a halo of the given 
mass are in gas and stars within the disks), is significantly larger 
over nearly the full range of Vrot- We speculate that the reason that 
the difference becomes less striking at higher \rot is that the Hi mass 
limit on which we select these galaxies is approaching the typical 
Hi mass for these more massive galaxies. 

In contrast, the star formation efficiency is remarkably similar 
for the Hi massive galaxies and typical galaxies, even as that 
efficiency appears to depend on Vrot- The rise in is interesting but 
difficult to interpret because it probably also reflects a change in 
the morphological mix of galaxies. It is nevertheless, particularly 
interesting that among the Hi massive galaxies, the efficiency fol¬ 
lows that of the galaxies despite variations of a factor of two 
in fc. This result suggests that star formation is strongly governed 
by the available gas mass and internal regulation, a result that is at 
least broadly reminiscent of results obtained in detailed studies of 



Figure 6. Deviations from the stellar mass Tully-Fisher relation, stellar, 
where the stellar TF is in effect the classic TF because the stellar masses are 
recovered directly from the luminosities. Panels and lines indicating means 
of the distribution are as in Figurej^ The deviation of the HI selected galax¬ 
ies is as strong as in Figure]^ even though these residuals are independent 
of HI mass. 

star formation rates in galaxies ( [Kennicutt & Evans|2012| l, rather 
than external factors. 


4.2 Environment 

Given our general understanding that massive galaxies tend to be 
early types, the physically interesting question is why the massive 
galaxies we have selected are such strong examples of late type 
galaxies. If the two primary drivers of galaxy evolutionary differ¬ 
ences among galaxies are mass and environment, and if mass is 
not the cause of the differences seen here, then environment has 
somehow played a dominant role for these galaxies. Confirming 
this conclusion is, however, difficult. By the nature of the selec¬ 
tion of the sample (isolated so that the Hi beam is clean and the 
result is a well-defined double horned profile), we have already se¬ 
lected against galaxies in dense environments. Upon examination, 
we confirm that none of these galaxies is in a highly overdense 
environment, such as an Abell cluster. Looking a bit further, we 
searched for differences in the environments of these galaxies rel¬ 
ative to galaxies in the 2MASS redshift catalog (2MRS) ( |Huchra| 
|et al.|201^ by comparing the correlation function of our galaxies 
to that of similar galaxies in the redshift catalog. We found no con¬ 
vincing evidence of a difference. 

We face two difficulties in exploring the role of environment 
further. First, environment means many different things. The local 
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Figure 7. The condensed baryon fraction, fc, which is defined as the frac¬ 
tion of all the cosmologically expected baryons in a halo that are present in 
the disk as either cold gas or stars, as a function of Vrot (upper panel) and the 
star formation efficiency, , which is defined as the fraction of condensed 
baryons that are in stars -i- Mgas ))• In this figure blue open circles 

represent the HI selected sample and black filled circles represent the S^G 
sample. Points only plotted if bin has at least 3 galaxies in it. Error bars 
represent error in the mean. 


environment, whether a galaxy is a parent or satellite, appears to 
play an important role in the observable characteristics of galax¬ 
ies ( [Berlind et al.|2005l [Tinker et aL|2012| l. On the other hand, the 
large scale (several Mpc) environment also appears to play a role 
( |G6mez et al.|2003} . Therefore, any one measure of environment 
will be a far from perfect diagnostic of the role of environment 
in evolution. Second, galaxies, particularly most massive galaxies, 
have inhabited many different environments over their lifetime. The 
environment we measure today, especially for those in denser than 
average environments, is not a complete measure of the environ¬ 
mental history of a galaxy ( |De Lucia et al.||2012| ). Both of these 
issues suggest that standard analyses may at best yield only sub¬ 
tle differences, even if environment does play a role in a galaxy’s 
evolution. 

To pursue this topic further, we explore a novel characteriza¬ 
tion of environment using the local shear field, as used for the La- 
niakea discovery { [Tully et al.|2014| ). A Wiener Filter reconstruction 
applied to the Cosmicflows-2 database ( |Tully et al.|2013| ) recovers 
the underlying 3D velocity field with an effective resolution of a 
few Megaparsecs ( [Courtois et al.|201^ . By taking spatial deriva¬ 
tives of the velocity field, we compute the dimensionless shear ten¬ 
sor. The tensor’s eigenvalues (Ti, A 2 , A 3 ) describe the strength (rel¬ 
ative to the Hubble expansion) of compression (positive values) or 


expansion (negative values) along the eigenvectors of the shear: ei, 
62 , and 63 . At each position in space, we calculate the three eigen¬ 
values of the velocity shear tensor. By ordering these eigenvalues 
from most positive to most negative, we set thresholds that capture 
four possibilities. Flows can be inward on all three axes, the con¬ 
dition of a cluster, inward on two axes and outward on the third, 
the condition of a filament, inward on one axis and outward on 
two, hence a sheet, or outward on all three axes, hence a void. We 
define boundaries around contiguous regions with the same shear 
properties and the contours outline the cosmic web as reconstructed 
by the V-web algorithm ( [Hoffman et al.|201^ . The statistical un¬ 
certainty of the shear eigenvectors and eigenvalues are controlled 
by means of an ensemble of constrained realizations. The spatial 
cosmography and reconstructed velocity field of the giant galaxies’ 
environment is shown in Figure 

We compute the number of galaxies in cells that are classified as 
knots, filaments, sheets, or nothing (voids or no signal) for two red- 
shift catalogs V8K ( [Courtois et al.|20r3] ) (the most complete red- 
shift catalog within 8,000 km/s) and 2MRS ( [Huchra et al.|2012| l. 
Some of the studied volume is without signal (for example the zone 
of avoidance). It is classified as a void in the V-web, but we do not 
have galaxies there. By confining the analysis to within 8000 km/s, 
our samples of Hi massive and fast rotators now contain 18 and 
17 galaxies, respectively, see Tables 2, 3 and 4. As seen in Tabl^ 
our Hi giant galaxies live preferentially in filaments: 10 out of 18 
(56%) while only 21% of all galaxies within the same volume lo¬ 
cated today in a filament. Filaments represent 10% of the volume, 
contain 21% of all galaxies, and 56% of the giants. Adopting a 
binomial distribution with the probability of a ’’true” output (fila¬ 
ment) defined to be 0.21, we find that the probability of having 10 
or more ’’true” outcomes is 0.00025. so the signal is significant at 
greater than 3cr confidence. In contrast, although our fast rotators 
also show a hint of a preference for filaments, there the result is 
only statistically significant at ~ Icr confidence. 


5 CONCLUSIONS 

We investigated the properties of galaxies that potentially challenge 
the current paradigm of galaxy evolution. Albeit rare, these massive 
disk galaxies run counter to the idea that mass alone determines a 
galaxy’s star formation rate or its morphology. Of course there are 
other previously known reasons to disfavor such an extreme view, 
but well-specified, carefully defined samples have the potential to 
highlight these discrepancies. We selected the 100 most Hi mas¬ 
sive galaxies and the 100 most rapidly rotating disk galaxies from 
Hi catalogs available for public use at the Extragalactic Distance 
Database that contains thousands of coherently measured Hi pro¬ 
files. 

To compare with the overall disk population, we use the bary- 
onic Tully-Fisher (BTF) relation as a benchmark. That relation 
holds broadly because a nearly constant fraction of the available 
baryons in each galaxy’s halo settle into the disk ( [Zaritsky et al.[ 
[2014b[ l. The classic Tully-Fisher (TF) relation then holds because 
for most galaxies a fixed fraction of those settled, or “condensed”, 
baryons turn into stars. Deviations from the BTF therefore indicate 
a variation in the fraction of condensed baryons, /c, which could 
signal interesting evolutionary deviations from the typical galaxy. 
Deviations from the TF relation could similarly signal deviations 
fc, or they could signal deviations in the star formation efficiency. 

Using the Hi data in combination with available 3.6//m pho¬ 
tometry from Spitzer observations, we place our two samples on 




























Giant disk galaxies 9 


Table 1. Environmental study using the Cosmic-V-web computed within a sphere of 8,000 km/s radius. Number and percentage (respectively) of galaxies 
located in structures dynamically classified as knots, filaments, sheets and voids. The most massive HI galaxies are preferentially located in filaments, with a 
frequency that is a factor of two larger than a random galaxy. 


Objects in knots in filaments in sheets in void cells 


HI massive 
Fast rotators 
2MRS galaxies 
V8k galaxies 
Vweb cells 


0 (0%) 
1 (6%) 
586 (2%) 
500 (2%) 
2459 (0.2%) 


10 (56%) 
5 (29%) 
8065 (21%) 
6408 (21%) 
112639(10%) 


4 (22%) 
8 (47%) 
18236 (47%) 
14468 (47%) 
451487 (42%) 


4 (22%) 
3 (18%) 
11483 (30%) 
9324 (30%) 
518007 (48%) 



-15000 -10000 -5000 O 5000 10OOO A 5000 

SGX (km/s) 



Figure 8. The environment is studied using the Velocity-Cosmic-Web (V- 
web). Regions gravitationally collapsing along one direction are sheets, 
along two orthogonal directions are filaments, along three orthogonal direc¬ 
tions are knots. Regions expanding in three directions are voids. Contours 
in grey, blue and red are showing regions classified respectively as sheets, 
filaments and knots. The top panel shows a 3D view. The bottom panel is 
a closer view of a slice of 1000 km/s width on Supergalactic Y axis. The 
galaxies in 2MASS redshift catalog (2MRS) are plotted with black dots. 
The cosmography is given with Laniakea on the left and Perseus-Pisces su¬ 
percluster on the right. The most massive HI galaxies (blue spheres) are 
located preferentially in filaments with a ~ 100% higher probability than 
random galaxies in two different redshift catalogs, V8K or 2MRS. For the 
fast rotators (red spheres) we cannot confidently identify a preferred envi¬ 
ronment. 


the BTF. We find both samples deviate systematically from the ex¬ 
isting BTF. The rapid rotators tend to lie below the BTF. We discuss 
a number of possible reasons for this deviation, including a sugges¬ 
tion that the scaling relation is not linear at the massive end. Such 
behavior is also seen in massive early-type galaxies ( [Bernardi et al.| 
|2011| ). The systematic deviation might also hint at bottom heavy 
IMFs in more massive galaxies, again as seen in massive early type 
galaxies ( |van Dokkum & Conroy|2010[|Cappellari et al.|2012| l. 

The Hi massive galaxies lie systematically above the BTF. If 
these deviations are interpreted as variations in /c, then we con¬ 
clude that there are galaxies in which fc approaches 1, or in other 
words the entire allotment of baryons for that halo have found their 
way to the disk. In spite of this, or whatever other process may be 
driving the deviation from the BTF, the star formation efficiency, 
when quantified as the fraction of disk baryons that are in stars, 
is the same in these systems as in more typical galaxies. We con¬ 
clude that the star formation efficiency is driven by internal self¬ 
regulation, while fc can be altered by external factors. 

Although mass is clearly not the driver of the morphological 
or star formation history differences that exist between these galax¬ 
ies and other similarly massive galaxies, we have not identified the 
driver. Our galaxies avoid the densest environments, but so do other 
galaxies that are not similarly extreme. When looking in detail for 
environmental differences, we found that the most massive HI de¬ 
tected galaxies are located preferentially in filaments. We present 
the first evidence of an environmental effect on galaxy evolution 
using a dynamical definition of a filament. 
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Table 2. The 20 largest Hl-mass galaxies 


PGC 

Name 

J2000 coordinates 

Vhel 
km s~^ 

km 5’“^ 

error on W^x 
km s~^ 

inclination 

degrees 

^max 

km s~^ 

MhiIMq 

10^^ solar mass units 

morph, type 

log(,d25) 

in log of 0.1 arcmin 

17625 

NGC1961 

J054204.8+692243 

3935 

658 

8 

47.0 

449.86 

4.46 

4.2 

1.65 

9560 

NGC0958 

J023042.8-025621 

5739 

570 

10 

78.1 

291.26 

3.77 

4.9 

1.40 

14030 

UGC02885 

J035302.4+353522 

5804 

556 

8 

62.4 

313.70 

3.35 

5.2 

1.63 

8066 

NGC0812 

J020651.6+443418 

5161 

435 

16 

71.9 

228.82 

3.24 

5.9 

1.44 

64458 

UGC11537 

J201838.1-000902 

5418 

499 

8 

69.7 

266.02 

3.09 

5.3 

1.18 

34836 

NGC3646 

J112143.1+201011 

4246 

508 

16 

63.9 

282.84 

2.96 

4.6 

1.49 

52328 

NGC5720 

J143833.3+504855 

7784 

409 

17 

52.1 

259.16 

2.80 

3.0 

1.28 

56891 

ESO136-016 

J160349.3-605840 

5426 

521 

12 

90.0 

260.50 

2.76 

5.2 

1.52 

24830 

UGC04625 

J085017.8+032951 

8476 

422 

18 

90.0 

211.00 

2.74 

6.1 

1.23 

62178 

NGC6674 

J183833.9+252230 

3429 

431 

16 

62.9 

242.08 

2.63 

3.0 

1.60 

16537 

UGC03218 

J050043.7+621439 

5228 

495 

15 

56.4 

297.15 

2.62 

3.1 

1.24 

59884 

IC4633 

J171346.9-773210 

2945 

318 

13 

62.4 

179.42 

2.57 

6.0 

1.43 

61791 

ES0395-002 

J182226.4-354040 

5607 

538 

13 

90.0 

269.00 

2.56 

3.4 

1.32 

19531 

NGC2280 

J064449.1-273819 

1897 

384 

7 

66.2 

209.84 

2.49 

5.9 

1.81 

7387 

NGC0753 

J015742.2-F355457 

4902 

314 

8 

52.5 

197.89 

2.41 

4.9 

1.14 

52361 

UGC09437 

J143910.6-F184247 

14401 

404 

16 

65.6 

221.81 

2.38 

3.3 

0.90 

54250 

NGC5833 

J151153.6-725134 

3031 

407 

14 

76.7 

209.11 

2.36 

4.2 

1.49 

9399 

NGC0931 

J022814.5-F311840 

4999 

422 

9 

81.3 

213.45 

2.30 

3.6 

1.39 

18089 

UGC03375 

J055525.3-F515438 

5791 

490 

11 

64.2 

272.13 

2.29 

5.2 

1.18 

28196 

NGC2998 

J094843.6-F440453 

4772 

374 

7 

61.8 

212.18 

2.27 

5.2 

1.39 
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Table 3. The 20 fastest rotators 


PGC 

Name 

J2000 coord. 

Vhel 
km s~^ 

Wmx 
km s~^ 

error on Wmx 
km s~^ 

inclination 

degrees 

^max 

km s~^ 

MhiIMq 

10^^ solar mass units 

morph, type 

log{d25) 

log of 0.1 arcmin 

17625 

NGC1961 

J054204.8+692243 

3935 

658 

8 

47.0 

449.86 

4.46 

4.2 

1.65 

42407 

NGC4594 

J123959.4-113723 

1095 

733 

17 

59.4 

425.80 

0.0430 

1.1 

1.93 

63575 

NGC6824 

J194340.7+560634 

3549 

544 

12 

45.8 

379.41 

0.148 

2.3 

1.28 

9359 

NGC0936 

J022737.4-010921 

1436 

545 

11 

50.5 

353.15 

0.0907 

-1.2 

1.65 

165398 

PGC165398 

J043157.1+592547 

4626 

489 

20 

45.5 

342.80 

0.440 

0.0 

0.94 

58596 

NGC6195 

J163632.6+390141 

9021 

517 

18 

50.2 

336.47 

0.221 

3.1 

1.16 

65375 

NGC6962 

J204719.1+001915 

4222 

475 

13 

45.0 

335.88 

1.12 

1.7 

1.43 

24685 

ESO563-021 

J084717.0-200208 

4582 

664 

7 

90.0 

332.00 

2.17 

4.3 

1.48 

26512 

NGC2841 

J092202.5+505837 

633 

592 

12 

65.2 

326.07 

0.236 

3.0 

1.84 

25161 

NGC2713 

J085720.5+025517 

3916 

629 

18 

77.7 

321.89 

0.503 

2.5 

1.53 

36706 

NGC3884 

J114612.2+202330 

6947 

462 

20 

45.9 

321.67 

0.660 

0.6 

1.28 

52665 

NGC5746 

J144456.0+015717 

1723 

634 

11 

90.0 

317.00 

0.456 

3.0 

1.86 

57173 

UGC10205 

J160640.2+300556 

6562 

541 

20 

59.0 

315.58 

0.725 

1.0 

1.16 

14030 

UGC02885 

J035302.4+353522 

5804 

556 

8 

62.4 

313.70 

3.35 

5.2 

1.63 

37617 

NGC3992 

J115735.9+532228 

1048 

459 

8 

47.4 

311.78 

0.218 

4.0 

1.91 

72233 

UGC12755 

J234349.7+282021 

8794 

503 

19 

54.7 

308.16 

0.648 

3.1 

1.07 

67966 

UGC11893 

J220406.7+355618 

5589 

606 

20 

80.2 

307.48 

0.924 

6.3 

0.81 

66880 

UGC11758 

J213057.6+135910 

8635 

612 

19 

90.0 

306.00 

0.942 

4.2 

1.12 

45947 

NGC5032 

J131327.0+274808 

6413 

524 

19 

59.5 

304.08 

0.212 

3.0 

1.29 

63286 

UGC11455 

J192956.3+720646 

5392 

608 

15 

90.0 

304.00 

1.29 

5.8 

1.39 
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Table 4. HI-21cm spectrum and Spitzer 3.6// images of the top 20 most heavy Hl-mass galaxies, in decreasing mass order. 


PGC - Name HI mass (logiMni) in solar units and Yhei PGC - Name 

Telescope / HI profile Spitzer 3.6 ju IRAC 400X400 pixels Telescope / HI profile 


HI mass (logiMni) in solar units and Yhei 
Spitzer 3.6 // IRAC 400X400 pixels 


PGC17625 



GET 

PGC14030 


PGC0014030-UGC028a5-ctf2009-GBT 



GET 


PGC64458 



GE300 


PGC52328 



GE300 

PGC24830 


PGC0D248J0-UGC04625-5hg2005-A0lf 


L lir\A.» A . r 


WVl J /~\ A \ 

fM tW yifi 






4.46/3935 itmi-' 


3.35/5804 



3.09/5418 kms-^ 


1 



2.80/7784 kms-' 



2.74/8476 kms-^ 


PGC9560 



Parkes 64m 


PGC8066 



GE300 


PGC34836 



GE140 


PGC56891 



Parkes 64 m 


PGC62178 



3.77/5739 Am i-' 


3.24/5161 Am 


♦ 



2.96/4246 Ami"* 



2.76/5426 Am 



2.63/3429 Am 



Arecibo 


GE140 
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Table 5. Table4 .. continued ...HI-21 cm spectrum and Spitzer 3.6// images of the top 20 most heavy Hl-mass galaxies, in decreasing mass order. 


PGC - Name HI mass (logiMni) in solar units and Yhei - PGC - Name 

Telescope / HI prohle Spitzer 3.6 // IRAC 400X400 pixels Telescope / HI profile 


HI mass (logiMni) in solar units and Vhei 
Spitzer 3.6 fi IRAC 400X400 pixels 


PGC16537 


2.62/5228 kms-^ 



GB300 



PGC59884 



Parkes 64m 




PGC19531 



Parkes 64m 


2.491 mi km s-' 



PGC7387 


PG[:0QQ73a7-UGCQ1437-shg2QQ5-AQIf 



l\A 


f...... ‘ 


v\J 


Arecibo 


2.41/4902 kms-^ 

- , ; ^--B—B 




b a. 


PGC52361 2.38 / 14401 km i”' 



PGC54250 


PGC0054250-NGC005S33-l<sk20Q4-PAKS 



Parkes 64m 


2.36/3031 kms-^ 



PGC9399 


2.3014999 km s-' 


PGC0aQ939g-UGCO1935-shg2pQ5-AQlf 









4600 4800 5000 5200 5400 


Arecibo 



PGC18089 



2.29/5791 kms-^ 




PGC28196 


2.21 lArlll kms-^ 


PGC0028T 

9B-UGC05250-sh(^2i 

DQ5-GB300 
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